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Abstract. It is well established numerically that spectral statistics of pseudo- 
integrable models differs considerably from the reference statistics of integrable and 
chaotic systems. In [PRL 93 (2004) 254102] statistical properties of a certain quantized 
pseudo-integrable map had been calculated analytically but only for a special sequence 
of matrix dimensions. The purpose of this paper is to obtain the spectral statistics of 
the same quantum map for all matrix dimensions. 
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1. Introduction 

The relation between spectral statistics of a quantum system and its classical counterpart 
is one of the main achievements of quantum chaos. It is established that at the scale 
of the mean level density the spectral statistics of classically integrable systems are 
described by the Poisson distribution [1] and the spectral statistics of classically chaotic 
systems are close to the statistics of eigenvalues of standard random matrix ensembles 
depending only on the underlying symmetry [2]. Though these statements had not 
mathematically been proved in the full generality and there exist noticeable exceptions 
(e.g. chaotic systems on constant negative curvature surfaces generated by arithmetic 
groups [3]), these conjectures are widely accepted for 'generic' quantum systems. 

Much less is known when a system is not classically integrable or completely chaotic. 
An important example which we have in mind is the case of pseudo-integrable systems 
(see e.g. [1]) represented by 2-dimensional polygonal billiards whose each angle is 
a rational multiple of tt. A typical classical trajectory in such models covers a 2- 
dimensional surface of a finite genus > 2 [5J. For comparison, a trajectory of a 2- 
dimensional integrable model belongs to a 2-dimensional torus of genus 1 and a typical 
trajectory of a 2-dimensional chaotic system covers ergodically the whole 3-dimensional 
surface of constant energy. 

Numerical calculations [6] clearly demonstrated that the spectral statistics of 
pseudo-integrable billiards is not universal and depends on billiard angles. It appears 
that spectral statistics of such models has two characteristic features: a level repulsion 
at small distances and an exponential decrease of the nearest-neighbour distribution 
at large separations. This type of statistics (called intermediate statistics) has been 
observed for the first time in the numerical simulations of the Anderson model at the 
point of metal-insulator transition [7j and later in this context has been profoundly 
investigated (see e.g. a recent review [8] and references therein). Unfortunately 
analytical progress in the investigation of intermediate statistics in pseudo-integrable 
systems is limited. In |9] the level compressibility of certain pseudo-integrable billiards 
was computed analytically, thus confirming the intermediate character of their spectral 
statistics. The main difficulty in analytical treatment of 2-dimensional pseudo-integrable 
billiards is the strong diffraction on billiard corners with angles ^ ir/n with integer n. 
Though there exists the exact Sommerfeld's solution for the scattering on such angles 
[To] , the diffraction coefficient is formally infinite along the optical boundaries so it 
is impossible to treat the multiple corner diffraction in perturbation series. A certain 
progress has been made in [H] where it was found that the multiple diffraction along 
periodic orbit channels in pseudo-integrable systems forces the wave functions to tend 
to zero on the channel boundaries thus forming superscarring states observed recently 
in micro- wave experiments [T2] . 

The difficulties in analytical solution of pseudo-integrable billiards lead to the 
necessity of investigation of simpler models with similar features. A promising example is 
the quantization of interval exchange maps which are known to be the correct description 
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of classical dynamics in pseudo-integrable systems [1], 

In [13] the following classical parabolic map had been quantized 

modi. (1) 



q + f{p + a) 

Here a is a constant and f{p) is a certain function (taken equal 2p in [13j). For rational 
a = m/q this map corresponds to the simplest interval exchange map of 2 intervals. 

After a straightforward generalization of the results of [13] the NxN unitary matrix 
associated with the above map has the form of the diagonal matrix e^** with real $a: 
multiplied by a constant unitary matrix fi^p {k,p = 0,1, N ~ 1) 

Mkp = e^*>fcp , (2) 



where 



N-l rr, ■ 



A^ 

r=0 



^ r{k — p + aN) 



1 _ 2-KiaN 

(3) 



jV(l - Q2-Ki{k-p+aN)/N-^ 



The matrix ^kp depends on a parameter a and the last equality are valid when aN is 
not an integer. (In the latter case the spectrum of ([2]) can be obtained analytically [2] 
and we always assume below that for rational a = m/q, mN ^ mod q.) 

Two main cases are of interest. The first corresponds to the ensemble of random 
matrices ([2]) where all A^ phases $fc are considered as independent random variables 
distributed uniformly between and 27r. We call such matrices non-symmetric ensemble. 

In the second case only a half of the phases are independent random variables 
uniformly distributed between and 27r. The other half is related to the first one by 
the symmetry 



This case will be called symmetric ensemble as in the dynamical interpretation the 
transformation k —>■ —k is the time-inversion symmetry. 

The eigenvalues A„ and the eigenvectors Uk{n), n = 1,2, . . . , N, are defined as usual 

Af-l 

Kuk{n) = y^^MkpUp{n) . (5) 

p=0 

Because matrix M^p is unitary {M M"^ = 1), its eigenvalues lie on the unit circle 

An = e'^" (6) 
and its eigenvectors can be chosen orthogonal 

N-l 

Uk{m)uk{n) = bmn (7) 

fc=0 

The statistical properties of matrix ([2]) depend crucially on the arithmetic of the 
parameter a. For irrational a. the map ([1]) is only a parabolic map and methods 
developed in [TS] and in the present paper can not be directly applied. Numerical 
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Figure 1. Open circles are nearest-neighbour distribution with a = \/5/4 and 
N = 801 for non-symmetric matrices (a) and symmetric matrices (b). Solid lines 
indicate the Wigner surmise ([8]). 



calculations suggest that for Diophantine a the spectral statistics of the matrix ([2]) is 
very close to the standard statistics of the Gaussian ensembles of random matrices |16j . 
Namely, non-symmetric matrices are described by the GUE statistics and symmetric 
matrices by the GOE statistics. For illustration in Fig. [1] we plot the nearest-neighbour 
distribution for matrices ([2]) with a = Vb/A and = 801. The solid lines indicates the 
Wigner surmise ([8]) which is known (see e.g. [16j) to be a good approximation for the 
GUE and the GOE distributions. 

Pgve{s) = ^s'e-'^'/- , pcoe(s) = Jse"-^'/^ . (8) 



Rational a = m/q with co-prime integers m and q correspond to an interval exchange 
map [13] and we shall consider only such values of a throughout the paper. 

This paper investigates the following question: What are the statistical properties 
of eigenvalues of the matrix ([21) for fixed rational a = m/q and large Nl It appears that 
to get a well defined limit in this case it is necessary to consider increasing sequences of 
matrix dimensions, A^, such that the product of the numerator of a times has a fixed 
residue modulo the denominator of a 

mN = k mod q . (9) 

In [15] it was demonstrated that for a special sequence of matrix dimensions with 

mN = ±1 mod q (10) 

eigenvalues of the matrix ([2]), ([3]) are described by the so-called semi-Poisson statistics 
which has been proposed in [17j as the simplest model of intermediate statistics. 

Let Xi < X2 < ■ ■ ■ < xk be an ordered sequence of real numbers (eigenvalues). The 
joint distribution for the semi-Poisson statistics is proportional to the product of the 
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nearest distances between these level times a confining potential V{x) 

P(xi, X2,...,xk)^1[ - J] e-^(^») (11) 

i i 

In the limit ii' ^ oo all correlation functions of the semi-Poisson statistics at the scale 
of the mean level density do not depend on V{x) and can be obtained analytically |17j . 
In particular, the probability that between 2 levels there exist exactly n — 1 levels has 
the form 

Pn{f3; s) = (^±l)!!^,n/^+-ie-(/^+i)^ . (12) 
T{np + n) 

The semi-Poisson statistics depends only on one parameter (3 which fixes the level 
repulsion at small distances so the nearest-neighbour distribution (i.e. Pn(/3; s) for n = 1) 
tends to zero as 

p{(3;s) = Aps^e-^^+^'>' (13) 

with (/5+l)^+Vr(/5 + l). 

From a mathematical point of view the semi-Poisson statistics can be considered as 
a stochastic process with independent increments (with gamma-distribution) forming a 
convolution semigroup 

{Pn*Pm){P;s) = Pn{P;y)Pm{P;s - y)dy = Pn+m{P;s) . (14) 



Jo 

According to [15] when the condition ( ITOi) is satisfied the spectral statistics of the matrix 
(j2]) tends for large to the semi-Poisson distribution with the following integer and 
half integer values of (3 depending on the denominator oi a = m/q and the symmetry 
of the map 

^ I g — 1 for non-symmetric ensemble 
I \q — 1 for symmetric ensemble 

To compare numerical calculations with theoretical predictions it is often more precise 
to consider instead of the nearest-neighbour distribution, p(s), (as in Fig. [1]) its integral 

N{s) = [ p{s')ds' (16) 



which gives the relative number of levels when the distance between the nearest- 
neighbour eigenvalues is less than s. This quantity is one of the main spectral correlation 
functions and throughout the paper we focus exclusively on it though other correlation 
functions can also be calculated and are of interest. 

To illustrate the convergence of the spectral statistics of the above matrices to the 
predicted values let us consider e.g. a = 1/2 with odd N. From (fTSl) and (fT3l) it follows 
that the limiting integrated nearest-neighbour distribution in the non-symmetric case is 
here the simplest semi-Poisson distribution with (3 = 1 

iV,p = 1 - (2s + l)e-2« . (17) 




Figure 2. Differences between the integrated nearest- neighbour distributions for the 
non-symmetric ensemble of matrices ([2]) with a = 1/2 for different odd N and the 
theoretical prediction for this case (HZ]). The different lines from bottom to top at 
small s correspond, respectively, to = 101, 201, 401, 801, and 1601. 



In Fig. [2] the difference between tlie integrated nearest-neiglibour distribution computed 
numerically and its tlieoretical prediction (fTTl) is plotted for different odd matrix 
dimensions. For this and other similar figures in the paper the number of realizations 
is chosen to be the minimum between 100 and 50000/iV. The figure shows that the 
agreement is quite good (less than 0.02) even for N of the order of a few hundreds. 

The purpose of this paper is to calculate the spectral statistics of the above 
quantized pseudo-integrable map in the general case (Q with k ^ 0,±1 mod q. The plan 
of the paper is the following. Sections [2] and [3] give extended details of the construction 
briefly discussed in [15]. The peculiarity of the problem under consideration is the 
existence of 2 rank-one deformations of the original matrix ([2]) with known eigenvalues 
and eigenvectors. These deformations are discussed in Section [21 In Section [3] it is 
demonstrated that these rank-one deformations lead to long-range correlations between 
the eigenvalues of the initial matrix ([21). To obtain a clear picture of these correlations 
it is convenient to use a special form of eigenvalue ordering (an analog of the unfolding) 
which is discussed in Section [H In Section [51 it is shown that these long-range correlations 
can effectively be taken into account by the construction of a kind of transfer operator. 
This operator is a finite dimensional matrix whose largest eigenvalue and corresponding 
eigenvectors permit to calculate all correlation functions. This is done explicitly for a 
few main examples in Section [61 Obtained analytical formulas agree well with numerical 
calculations. The summary of the results is present in Section [7l Certain technical 



details are given in Appendix A 
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As was shown in [15] the important property of the matrix (Ej), ([3]) is the possibility to 
rewrite it in the following form 



f, = e^-fft ^- - ^ = m + - - ( 



k 



where a new matrix Nkp is 

Nkp = Mfcpe2"'('=-P+"^)/^ . (19) 
Eigenvalues A'^ and eigenvectors ipk{f^) of the matrix Nkp 

N-l 

KMn) = Y.^>'pMn) (20) 

can easily be expressed through the eigenvalues and the eigenvectors of the original 
matrix M^p 

Mn)=e'''''/''uk{n), A^ = e^^A. . (21) 

But from f|T8l) it follows that matrix Nkp is a rank-one deformation of matrix Mkp so 
it is possible to construct its eigenvalues and corresponding eigenvectors in a different 
way (see e.g. [IS])- 

Write a formal expansion of an eigenvector of matrix Nkp as a series of the complete 
set of eigenvectors of the matrix Mkp 

N 

Mn) = e'-''/''uk{n) = ^ c^{n)uk{m) . (22) 

m=l 

From f|T8l) one gets 

K Cm{n)ukim) = ^ Cmin)AmUkim) — e'** ^ 0^(72) ^ Up{m) . (23) 

m=l m=l m=l p=0 

Introducing the notations 

N-l N 

Am = ^Mfc(m) , g{n) = ^Cm{n)Am (24) 

k=0 m=l 

and using the orthogonahty of eigenvectors Uk{n) ([7j) one obtains 

Cm{n) = 9{n) " (25) 

where 

5^ = 5Ze'*^nfc(m). (26) 

k=0 

Multiplying the both sides of fl25l) by and summing from 1 to iV one concludes that 
every eigenvalues A^ of the matrix Nkp obey the equation 

1 „27riaAf ^ /I R 

A - A' • ^^^^ 



iV ^ A^ - A 



m=l 



Spectral statistics of a pseudo-integrable map: the general case 



8 



From and it follows that 

ArnB*^ = (28) 
and Eq. (127|) takes the final form 

m " 

For rank-one deformations of a real symmetric matrix all terms in the corresponding 
equation would be real and one easily comes to the well known conclusion that 
eigenvalues of a rank-one deformation of a real symmetric matrix are in-between 
eigenvalues of the unperturbed ones (cf. [IBD- Iii our case the both matrices, Mkp 
and Nkp are unitary and their eigenvalues lie on the unit circle: = e^^", = e'^". 
So arguments require straightforward modifications. 
From (129|) one gets 

\ — ^ |-SjTlP 

m 

which can be rewritten as follows 

N{cot TiaN -i) = J2\Bj^ (cot ^" ~ - i) . (31) 

Due to the completeness of Uk{rn) one has 

N 

^Up{n)ukin) = Spk . (32) 

n=l 

and, consequently, 

N 

Y,\Bn.\' = N. (33) 

771=1 

Therefore, the imaginary part of fl3T]) is identically zero and new phases 6'^ have to be 
determined from a real equation 

F{e'J = N cot IT aN (34) 

where F{9) is defined by: 

F{e) = J2\B,n\'cot-^. (35) 

m=l 

In the interval [0, 27r) -F(^) has poles at 6 = 6m (assuming that all B^, 7^ 0) and it is 
monotone between them (cf. Fig. [H]). Let the eigenphases 6m be ordered on the unit 
circle < 6*1 < 6^2 < . . . < ^at < 27r. Then between two near-by eigenvalues 6m and 6m+i 
(mod 27r) there exists one and only one new eigenvalue Here the new eigenvalues are 
not necessarily ordered. But according to (12T|) all eigenphases of the matrix N^p have 
the form 6'^^ = 0m + 2vra; (mod 27r). Therefore we prove the following lemma (cf. Fig. H]). 




Figure 3. Schematic plot of Eq. ((Ml) (solid black line). F{e) is defined by p5|) . 
Vertical dashed lines represent the pole positions. Horizontal straight line indicates 
the value of the right-hand sine of Eq. (p4)l . The abscissa of its intersections with F{d) 
(indicated by black circles) give the solutions of that equation. 




Figure 4. Illustration of Lemma 1 . Black circles denote the position of 7 eigenphases 
for a — 1/5 and black lines are their radius- vectors. Black squares indicate the position 
of eigenphases after the rotation by angle 2n/5 and dashed red lines are radius- vectors 
of the rotated eigenphases. The rotated points are indicated by the same number but 
with sign '. 



Lemma 1. The eigenvalues of the unitary matrix M^p defined in ^ 
and 1^ are such that after the rotation by 27ra in-between of any pairs 
of nearest eigenvalues there exist one and only one rotated eigenvalue. 

Multiplying ( ITSll by exp(— 27ri(fc — p + aN)/N) one gets another relation 

1 _ p27riaAf 
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with a new matrix N^p 

Nkp = Mfcpe-'"'(^-^'+"^)/^ (37) 
whose eigenvalues and eigenvectors are 

V^fc(n)=e-2-'^/V(r2), = e-2™A„ . (38) 

Repeating the above calculations but for the matrix N^p one finds that its eigenvalues 
have to be determined from the equation 



1 _ p2vriaAf ^ A If) 12 
1 e V V i\.'m.\-D^\ 



N ^ A^ - A' 



(39) 



where 

N-l 

5™ = $^^fc(m)e^*'=-2-'=/^ (40) 

As it has exactly the same form as ( l29l) and A^ = e~^'^'°A„ one comes to the lemma. 

Lemma 1'. The eigenvalues of M^p are such that after the rotation by —2iTa be- 
tween two nearest eigenvalues of M^p there exists one and only one rotated eigenvalue. 

These lemmas prove the existence of long-range correlations between eigenvalues of 
the matrix Mkp which are merely a consequence of the fact that rank-one deforma- 
tions ( JTSl) and ( |36l) of the original matrix ([2]) have eigenvalues easily expressible through 
eigenvalues of the original matrix. 

A few other consequences of this property is worth to mention. As all solutions 
of fl2^ have the form A^ = A^e^'^'" with n = 1, . . . ,N the numerators of this equation 



can be found explicitly. From Appendix A it follows that 



TV 

rr(A„-A„e2™) 

Am\Bm\'' = (41) 

which can be rewritten in the real form as follows 

^ simraN ^ sin(|(g^ - g„ - 27ra)) 

' Nsinna ^} smU (Om - On)) ' ^ ' 

Similarly, from (15^ one concludes that 

N 

n(A™-A.e--) 
-A^lfi^P = . (43) 



N 



Yl {Am - A. 



So 



~ ^ sinvraA^ _ sin(|(6'n, -9n + 2na)) 
' ATsinvra ~ ^} sin(i(^^ - 0„)) ' ^ ^ 
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n 



k 



k 




Figure 5. (a) Eigenphases close to the boundaries of a 27ra sector, (b) Division 
of the unit circle into sectors of angle 2-Kra/q. Black circles indicate the position of 
eigenphases in two near-by sectors. Dashed lines show the positions which will occupy 
eigenvalues from one sector after the rotation by 2-Km/ q. 



3. Long-range correlations 

In the preceding Sections it has been proved that eigenphases of matrix ([5]) have a 
special type of long-range correlations. Namely, when one rotates all eigenvalues of 
this matrix by ±27ra and superimposes the rotated eigenphases with non-rotated ones 
then in-between two nearest eigenvalues of the original matrix there will be one and 
only one rotated eigenphase. In this Section we investigate certain consequences of such 
correlations in more details. 

Let us put all eigenvalues of unitary matrix ([2]) on the unit circle and consider a 
sector of angle 27ra which contains n eigenvalues (see Fig. [5^). The sector boundaries 
divide the unit circle in-between certain eigenphases. Denote the angular distance from 
the sector boundaries to the nearest eigenphases in the clockwise and counterclockwise 
directions by y^, yk+i and x^, x^+i respectively (cf. Fig. [5^). After the rotation by 27ra 
only 2 possibilities are possible: either the point Xk or the point yt will fall inside the 
points Xfc+i and Hk+i- In the first case one has < x^+i and yk+i < yk- In the second 
case the inequalities are reversed: Xk > Xk+i and yk+i > yk- Therefore in the all cases 
the following inequality is fulfilled 



This inequality is valid for all a and A^. From now on we shall consider only rational a 



with co-prime integers m and q. 

As above divide the unit circle into q radial sectors of angle 27cm/q. When m = 1 
these sectors are disjoint but for m > 1 they will overlap. Denote the number of 




(45) 



m 

a = — 



(46) 
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eigenphases in the k^^ sector by Uk (see Fig. [5b). 

After the rotation by 2'Km/q the eigenphases from the sector will move into 
the (fc + 1)^^ sector. These rotated points will divide this sector into + 1 intervals. 
According to the lemma 1 the eigenphases in the (fc + 1)^'^ sector have to be intertwined 
with the rotated eigenvalues. Therefore, all intervals except the first and the last have 
to be occupied. The first will be occupied if Xk > Xk+i and the last will be occupied if 
Vk+i > yk+2- AH these requirements can be rewritten as the following recurrence relation 

Hk+i = Uk - 1 + Q{xk - Xk+i) + Q{yk+i - yk+2) ■ (47) 

Here as in Fig. [5^ Xk and yk are distances from the boundary of the k^^ sector to the 
two closest eigenphases to it and 9(x) is the Heaviside function: 9(a;) = 1 if a; > and 
9(x) = if X < 0. From ( H5l) it follows that the difference yk+i — yk+2 is of the opposite 
sign than the the difference Xk+i — Xk+2 and the last relation takes the form 

Uk+i = Uk - 1 + Q{xk - Xk+i) + Q{xk+2 - Xk+i) ■ (48) 

Now choose the beginning of the first sector at the position of an eigenphase (i.e. impose 
yi = 0). Then from fHSj) it follows that X2 < Xi and Xq < Xi. Direct applications of (jUj) 
give 

712 = rti - 1 + 0(a^i - X2) + 6(a:;3 - X2) = ni + 0(3:3 - X2) (49) 

because Xi > X2, 

ns = n2 - 1 + Q{x2 - X3) + 6(^4 - X3) 

= ni + Q{x3 - X2) - 1 + 6(x2 - X3) + Q{x4 - X3) = ni + Q{x4 - X3) (50) 

because 0(x) + 0(— x) = 1 and so on. In this manner one concludes that for 
j=2,...,g-l 



n 



'0 



rii + 0(xj+i - Xj) (51) 



and rig = rii + 1 because, as was noted above, Xq < xi. 

As the sum over all q sectors cover the unit circle m times, the sum over all rik 



equals mN: Uk = mN. Therefore 



fe=i 

q-l 

mN = gm + 1 + ^ <d{xj+i - Xj) . (52) 

When mN = 1 mod q, all 0-functions in the right-hand side of this expression are 
forced to be zero which leads to the conclusion that in this case 

Xq < Xq-l < . . . < X2 < Xl . (53) 

When mN = — 1 mod q, all B-functions have to be equal to 1 and the inequalities are 
reversed 



X2 < X3 < . . . < Xq-l < Xq < Xl 



(54) 
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Figure 6. (a) Small black circles: schematic view of eigenvalues of matrix ^ for 
a — 1/5 and = 14 (N = — 1 mod 5). Numbers from 1 to 14 indicate the consecutive 
eigenphases. Eigenphases rotated by 2tt/5 are denoted by red squares. For clarity they 
are situated on a smaller dashed circle, (b) The same configuration of eigenvalues but 
unfolded on 5 horizontal lines representing 5 sectors, xi, . . . are the distances from 
the beginning of the sectors to the closest eigenphase as in Fig. [5^. Thick red lines 
demonstrate relative eigenvalue positions. 

Inequalities (153|) and ( l54l) manifest the existence of an exceptionally strong long-range 
correlations when mN = ±1 mod q. For usual matrix ensembles correlations between 
eigenvalues separated by a large distance tend to zero. But in our case eigenphases at 
distances 27rmfc/g with k = 1, . . . , (g — 1) are not independent but restricted by the 
above inequahties. In [15] only these special cases had been considered. 

4. Geometrical unfolding 

To visualize better the restrictions implied by inequalities (H5l) and recurrence relations 
f HHj) let us split the unit circle into q angular sectors of angle 2nm/q as above and denote 
the positions of eigenvalues inside each sector on q parallel lines (see Fig. [6]). 

According to Lemmas 1 and 1', between two closest eigenphases of the matrix 
([2]) there is one and only one eigenphase rotated by 27ra = 27rm/g. In the unfolded 
description it is manifested by the condition that points at each line have to be in- 
between two close-by points on the lower line. The relative positions of eigenvalues 
strongly depend on distances Xk from the beginning of each sector to the eigenphase 
closest to it (see (l52l)). In the unfolded representation (as in Fig. EK) Xk are the 
distances along the horizontal lines to points closest to a vertical line which represents 
the boundary of a sector with angle 27rm/g. 

Let start from the lower left point and draw a horizontal line till there is at the 
vertical a point situated at the first line above it. Then draw the vertical line till it 
touches that point. Now continue drawing an horizontal line till there is at the vertical 
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Figure 7. (a) Small black circles: schematic view of eigenvalues of matrix ^ for 
a ~ 1/5 and = 18 (i.e. N = —2 mod 5). Thick red lines connect eigenphases 
as indicated in the text, (b) The same configuration of eigenvalues but connected by 
non-decreasing staircase lines. 



a point at the upper line closest to the boundary of the given sector and so on. This 
line will go right if Xk+i > Xk and left if Xk+i < Xk- Finally, points will be connected 
by step-wise lines as in Fig. |6]d. Notice that according to our convention point 1 does 
not belong to the first line but to the last one. The shape of these lines are determined 
by the inequalities between all Xk- In Fig. [6] the case of = — 1 mod 5 is indicated. 
According to (15^ Xk+i > Xk for k = 1, ... ,4 which explains the staircase form of these 
lines. They all start at points along the first horizontal line, go up and to the right, and 
finally finish at the last horizontal line but with the shift by 1 unit. It is clear that such 
lines cannot cross each other. 

For other matrix dimensions these lines will have a different shape. Consider first 
as an example the case a = 1/5 and = 3 mod 5. From (!52|) it follows that in this 
case only 2 of 3 G-functions have to be equal to 1. Every time one of the 6-function is 
zero, the above lines turns to the left (see Fig. [7]) so the shape of the step-wise curve 
differs from the one of Fig. Eb- Consider the horizontal line when it first turns left. 
Instead of the left turn let us continue to the right till we touch an other point on this 
line. As between these two points there exists only one point at the lower and higher 
lines there is no contradiction with our line construction (cf. Fig. [7j). Finally we come 
to the conclusion that for these values of a and A^ the eigenphases have to be connected 
by non-intersecting lines which go only up and to the right and whose initial and final 
points are shifted by 2 units. 

These arguments can be generalized for all values of matrix dimensions and we get 
the following lemma. 

Lemma 3: For a = m/q and mN = k mod q with k = 1, ... ,g — 1 mutual posi- 
tions of eigenphases of matrix ^ can be described as follows. Fix q horizontal lines, 
put arbitrary points at the lowest line, and notice the vertical images of these points 
along the last line. Draw staircase non-intersecting lines going only up and to the right 
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with the condition that they start at the lower line and end at last line but with the shift 
by k units. Points at horizontal lines are situated at the corners of the constructed lines. 

When k > q/2 one may simplify the construction by using non-intersecting step- 
wise lines going up and to the left with the shift by g — A;. It implies that properties of 
the cases mN = k mod q and mN = —k mod q are the same. 

5. Transfer operator 

The above unfolding gives not only a clear picture of mutual positions of eigenphases 
but also serves as the basis of the explicit calculation of the spectral statistics for the 
problem under consideration. The calculations are based on the following conjecture 
proposed in [15]. 

Conjecture: For a = m/q the eigenvalues of the q^^ power of the original matrix 

for all N ^ mod q have universal spectral statistics independent on q and N 
but different for different symmetry classes of random phases For non-symmetric 
ensemble the statistics of coincides with the Poisson statistics and for symmetric 
ensemble 0) it is described by the semi-Poisson statistics with (3 = —1/2 called in IW^ 
the super-Poisson statistics. 

The main physical argument in favor of this conjecture in [1^ was the fact that 
for rational a = m/q the g*^ power of the classical map ([1]) corresponds to a classically 
integrable map, and according to the usual wisdom P integrable models have to be 
described by the Poisson statistics. Extensive numerical calculations agree very well 
with this conjecture. But it seems that to prove it rigorously one has to develop new 
methods which are at present under investigation [19]. 

After unfolding, when the points from all the sectors separated by 2Txm/q are taken 
together, they can be considered as the eigenphases of the g*^ power of the original 
matrix (with evident rescaling). Assuming the validity of the conjecture it follows that 
for all ^ mod q these points constitute the semi-Poisson ensemble with /? = for 
non-symmetric matrices and (3 = —1/2 for symmetric ones. 

In particular, the probability that between two eigenvalues of M"? separated by x 
there exist exactly r eigenvalues is the following 



The results of the preceding Section can be reformulated such that the joint probability 
of the close-by levels integrated over all the possible positions of other levels is the same 
as the probability of non-intersecting staircase lines which start from the initial levels 
and which finish after q steps (where q is the denominator of a) with the shift of k units 
(where k is the residue of mN modulo q). According to the conjecture the distribution 
of unfolded points are known which permits the calculation of spectral statistics of the 
original matrix ([2]). 




for non — symmetric ensemble 



for symmetric ensemble 



(55) 
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The usual method of computing the probabihty of non-intersecting paths in a 
Markoff process is the determinantal representation [20]. We found that for practical 
reasons it is more convenient to use the transfer operator method. In this method one 
first unfolds the spectrum as indicated in Fig. [71 Then one considers a vertical strip 
bounded by vertical lines emanating from any two near-by levels of original matrix, say 
X2 and X3 in Fig. [71 separated by the distance x = Xs — X2- Now different horizontal lines 
can enter and can leave the strip. When all in-coming and out-coming horizontal lines 
are fixed, it is obvious that the configuration inside the strip is not affected by outside 
points. Therefore, it is possible to integrate over all configurations of points inside the 
strip with prescribed ordering. Denoting the initial and final lines by indexes i and j the 
result of integration constitutes the matrix element Tji{x) of the transfer matrix T{x). 

For rational a = m/q and mN = k mod q with 1 < k < q — 1 each initial (and 
final) state is determined by fixing I = q — k — 1 horizontal lines 

l<ii<i2<...<ii<q-2 (56) 

from the total number of lines equal q — 2. It means that the dimension of the transfer 
operator is 

t = C\_, (57) 

and it is convenient to label the set of / integers obeying ( l56l) in e.g. lexicographical 
order. 

For clarity, let us consider the case a = 1/5 and = 3 mod 5 in detail (/ = 1). As 
g = 5 there is one possible horizontal line which may go through the vertical strip (cf. 
Fig. [7]). Therefore the transfer matrix is 3 x 3 matrix labeled by the number of these 
lines. In Fig. [HI all possible configurations for this case are presented. 

In general, iii = {ii,i2, ■ ■ ■ and j = (ji, 72, • • • ,ji) are multi-indexes of an initial 
and a final states, the total number of points, r, inside the considered vertical strip is 
determined by the expression 

i-i 

r = ji + ^[j.+i - is] + q - I - k = \j\ - \i\ + q - I (58) 
where \j\ = ji + . . . + ji and \i\ = ii + . . . + ii. In general, rmin < < '^max with 

''"min = 2 , Tjjiax = k(^q — k) . (59) 

If at least one term in the square brackets in fl58l) is negative, the configuration is 
impossible and the corresponding matrix element equals zero. Otherwise, the integration 
over all intermediate configurations compatible with the imposed inequalities gives the 
value of the transfer matrix elements. 

The calculation of this probability is straightforward. According to the above 
conjecture, the probability that between 2 eigenvalues of separated by x there exist 
r ordered eigenvalues i/s such that 

0<yi<y2<...<yr<x (60) 
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Figure 8. Structure of the transfer operator for a = 1/5 and = —2 mod 5. 



is given by fl55l) . Therefore, the transfer matrix element Tji{x) is the product of 2 factors 

Tji{x) = UjiPrix) (61) 

where r is the integer determined by (!58|) . Pr(a;) is the same as in (l55ll . and nji is the 
number of configurations of r points which fulfilled all inequalities comparable with the 
fixed initial and final states. Interchanging the initial and final states and counting 
horizontal lines from the top one gets that the transfer operator matrix elements obey 
the following symmetry 



Tji{x) = TiTjT^x) 



(62) 



where if i = (ii, . . . , i;), = (g - 1 - i;, . . . , g - 1 - ii). 

For example, for a = 1/5 and = —2 mod 5 the T12 element includes 5 points 
(cf. Fig. [8]), 3 points, a, 6, c belong to the upper curve, and 2 points, A and B, belong 
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to the lower curve. From the mutual positions of these points it follows that the Ti2{x) 
matrix element equals the probability that the following inequalities are fulfilled 

0<a<b<c<x, < A< B < X, a < A, b< B. (63) 

By inspection one finds that there exist exactly 5 possible ordered sequences compatible 
with inequalities (1631) . namely 

a b c A B, a b A B c, a A b B c, a b A c B, a A b c B. 



Therefore Ti2{x) = 5p5(x). Following the symmetry (I62I) we also have 
^23(3;) = Ti2(x) = 5p5(x). 

By construction the joint probability of the eigenvalues is equal to the product of 
the transfer matrices over all near-by points. As we are interested in the limit of large 
number of eigenvalues, the exact behavior near the boundaries are not important and 
one can take simply the trace of the whole product. Finally the joint probability of 
near-by levels of the original matrix < Xi < 0:2 < . . . < xk < L integrated over all 
possible configuration of levels on other sectors takes the form 

Pl{xi, . . . ,x/f) ~ Tr [T{xk - xk-i) ■ ■■■■ T{x2 - Xi))]5{xi+ X2 + . . . + xk - L) .(64) 

The further steps are standard (see e.g. |I3)- One has first to calculate the Laplace 
transform of the transfer matrix 



00 

kx 



T{k) = / T(a;)e-'=^da; . (65) 







Let A(A;) be the largest eigenvalue of T{k), and w{k), v{k) be respectively the right and 
left eigenvectors of this matrix corresponding to \{k) 

f{k)w{k) = X{k)w{k) , 'v{kyf{k) = X{kyv{k) . (66) 

In the limit K ^ 00 the dominant contribution comes [T7] from a vicinity of the saddle 
point, kgp = h, defined from the condition of the fixed mean level density L/K which 
we normalize to 1 

X'(h) , , 

^+1 = 0, (67) 

Then the nearest-neighbour distribution is determined by the formula 
Condition fl67|) is equivalent to the standard normalization of p{s) 







p{s)ds = 1, / sp{s)ds = 1. (69) 

From (1551) it follows that 

^ I (fc + l)^^'"+^\ for non — symmetric ensemble , , 

■'^ ' {2k + l)"(^'+i)/2, for symmetric ensemble 
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As r is determined by (1551) . the dependence of eigenvalues and eigenvectors (1^ on k is 
easy to find 

Wi = Wi{k + 1)-^'^ , Vi = Vi{k + lf\, X{k) = ~X{k + (71) 

for non-symmetric matrices, and 

w,{k) = iI;,(2A; + l)-l*l/2 , v,{k) = v,{2k + iy'^/^ , X{k) = X{2k + l)-''/\72) 

for symmetric matrices. Here i denotes the multi-index (ii, . . . ,ii), \i\ = ii + . . . + ii, 
and all tilded quantities are independent on k. 

From these relations one finds that the saddle point h obeying ( l67l) is 



q — 1 for non — symmetric ensemble 
(g — l)/2 for symmetric ensemble 



(73) 



Using ( l59l) we conclude that the nearest-neighbour distribution (l68i) for a = m/q and 
mN = k mod q equals the following finite sums 

k{q-k) 

p{s) = «ns"e-^^ (74) 

n=2 

for non-symmetric matrices and 

{k{q-k)-l) 

p{s)= Yl «n/2^"/'e-'^^/2 (75) 

n=l 

for symmetric ones. 

The nearest-neighbour distribution for all considered cases (with k ^ 0, ±1 mod q) 
manifests level repulsion at small s 

, , 1-5^ for non — symmetric ensemble , , 

I s^/^ for symmetric ensemble 

and has the exponential decrease at large s as it should be for intermediate statistics. 

Other correlation functions can also be written explicitly through the same 
quantities [17]. In particular, the two-point correlation form factor has the following 
form 

K{t) = 1 + 2Re ^(27rir) (77) 

where 

^ 'w{h)L{t^h){l-L{t^h))-^v{h) 
^ ^ ^w{h)v{h) ^ ' 

and the matrix L{s) = T{s)/X{h). 

One can check that for all the level compressibility K{0) = l/q for non-symmetric 
matrices and K{0) = 2/g for symmetric ones. 

Numerically it was estabhshed pi] but not yet proved analytically that eigenvectors 
of the considered ensembles of random matrices have fractal properties independent on 
the residue k ^ mod q. 
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The simplest new case corresponds to a; = 1/5 and = ±2 mod 5. Considering all 
configurations in Fig. [8], one gets that in this case the transfer matrix has the following 
form 

3p4{x) 5^5 (x) 5pq{x) ^ 
T{x) = Spsix) 5pi{x) 5p5{x) . (79) 

y 2^2(3;) 3^3(0;) 3^4(3;) J 

Performing the calculations discussed in the precedent Section we find that when 
a = 1/5 and = ±2 mod 5 the nearest-neighbour distribution for non-symmetric 
matrices is 

p{s) = (a2S^ + ags^ + 045^ + a^s^ + aGS^)e~^' (80) 

where coefficients a„ are the following: 02 = 625/2 — 275v^/2 ^ 5.041, 03 = 
3125/2-1375^/5/2 25.203, 04 = 71875/48 + 33125V^/48 ^ 45.724, ag = -15625/3 + 
9375^/5/4 ^ 32.451, ae = 1015625/288 - 453125\/5/288 ^ 8.357. 

In a similar manner one finds that for symmetric matrices under the same conditions 
p{s) is given by the following expression 

pis) = (ai/2S^/^ + ais + a3/2S^/^ + ags^ + 05/2^^/^)6"^" (81) 

with ai/2 ^ .3597, ^ 1.5122, 03/2 ^ 2.6105, 02 ^ 1.9471, 05/2 ^ .5725. 

For a = 1/7 and N = ±2 mod 7 the transfer operator is represented by the 5x5 
matrix: 

28p8(a;) 
28p7(x) 
19p6(a;) 
10^5(3;) 
4^4(0;) 

Computing its largest eigenvalue and using ( l68l) one finds that the nearest-neighbour 
distribution in this case for non-symmetric matrices has the form 



T{x) 



( ^Pq{x) 14^7 (x) 

5]95(x) 14p6(a;) 

4^4(0;) 10^5(0;) 

3^3(0;) 6^4(3;) 

\ 2^2(3;) 3^3(0;) 



42p9(x) 42pio(x) \ 

42p8(a;) 42^9 (x) 

28p7(a;) 28p8(a;) 

lAp^^x) 14^7 (x) 

5p5(x) hpfs{x) J 



(82) 



p{s) = {a2S^ + a^s^ + a4S^ + a^s^ + a^s^ 
+ 07^^ + a^s^ + ags^ + aioS^°)e~'^'' 



(83) 



where coefficients a„ are: 02 — 3.4998, 03 ~ 24.4986, 04 ~ 82.4309, 05 
176.8723, as ~ 251.6396, ~ 229.5488, ag ~ 130.8981, ag ~ 43.7932, aio ~ 6.8214. 
For symmetric matrices for the same a and N 

(ai/2v^ + dis + a3/2S^/^ + a2S^ + a^j^s^l"^ 



pis) 



+ ass^ + a7/2S^/^ + a4S^ + ag/2S^/^)e-^'^^ 



U) 



with ai/2 ~ .1508, ai ^ .7500, 03/2 ~ 2.0293, 02 
5.0193, 07/2 ~ 3.1567, ~ 1.2312, 09/2 ~ 0.2350. 



3.8675, 05/2 ~ 5.3099, 
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Figure 9. One of transfer matrix elements for a = 1/7 and N = —3 mod 7. 



For a = 1/7 and = ±3 mod 7 there exist two possible entering lines and two 
leaving lines (cf. Fig. [9l) The dimension of the transfer matrix in this case is C| = 10. 
Its explicit form is the following 



r(x) 



/ 10P6 


35p7 


70p8 


84p9 


56p8 


168p9 


252pio 


210pio 


462pii 


462pi2 


\ 


10p5 


35p6 


70p7 


84p8 


56p7 


168p8 


252p9 


210p9 


462pio 


462pii 




Qpi 


20p5 


401^6 


A9p7 


30p6 


91p7 


140p8 


112p8 


252p9 


252pio 






8P4 


15p5 


WPQ 


10p5 


30p6 


49p7 


35p7 


84p8 


84p9 




4p4 


15p5 


30p6 


35p7 


26p6 


77p7 


112p8 


98p8 


210p9 


210pio 




3P3 


12p4 


25p5 


30p6 


20p5 


61p6 


91p7 


77p7 


168p8 


168p9 




2P2 


6P3 


12^J4 


15p5 


8P4 


25p5 


40p6 


30p6 


70p7 


70p8 







3p3 


8P4 


WP5 


6p4 


20p5 


30p6 


26p6 


56p7 


56p8 









6P3 


8P4 


3P3 


12p4 


20p5 


15p5 


35p6 


35p7 




V 





2P2 


3P3 





3P3 


6p4 


4p4 


WP5 


Wp6 


I 



15) 



Finally one obtains that for a = 1/7 and = ±3 mod 7 the nearest-neighbour 
distribution for non-symmetric ensemble is 

p{s) = (a2S^ -|- a^s^ + a^s'^ + a^s^ + a^s^ + ajs''' 



+ ass^ + ags^ + aios^^ 



aiis 



,11 



«i2S^^)e" 



■7s 



4.056, ag ~ 28.3898, ^ 91.6591, ~ 177.9134, ag ~ 227.8782, 



where 02 

200.0096, 08 ~ 121.6091, ag ~ 50.5880, aio ~ 13.778, an ~ 
For symmetric matrices under the same conditions 

p{s 



(86) 
07 ~ 



2.2159, ai2 ~ .1596. 



7s/2 



^7) 



with ai/2 ~ .1747, ai ~ .8691, 03/2 ~ 2.2565, 02 ~ 3.8902, 05/2 ~ 4.8085, ~ 
4.3734, 07/2 ~ 2.9327, ~ 1.4222, 09/2 ~ .4747, 05 ~ .0979, au/2 ~ .0094. 

In Fig. [in] and [11] the calculated nearest neighbour distributions are plotted for 
a = 1/5 and a = 1/7 with all possible residues of ^ modulo 1/a. As expected, 
the case A^ = ±1 mod 1/a differs considerably from other cases. When the residue, k, 
increases the nearest-neighbour distribution more and more resembles to the nearest- 
neighbour distribution of the standard Gaussian ensembles of random matrices. For 
example, for a = 1/7 the results with A^ = ±2 and A^ = ±3 are difficult to distinguish 
from the Wigner surmise ([8]) for GUE (for non-symmetric matrices) and for GOE (for 
symmetric ones) 
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a) b) 

Figure 10. Nearest-neighbour distribution for a — 1/5 for (a) the non-symmetric 
ensemble and (b) the symmetric one. Dashed lines correspond to N = ±1 mod 5 given 
by (fT^ with f3 — 4 in (a) and P — 3/2 in (b). Solid lines indicate the results for 
N = ±2 mod 5 given by ^ in (a) and by ^ in (b). 




S 



a) b) 

Figure 11. The same as in Fig. [10] but for a = 1/7. Dashed black lines correspond 
to iV = ±1 mod 7 given by (jl3p with /3 = 6 in non-symmetric matrices and /? = 3/2 
for symmetric matrices. The red dotted dashed lines indicate the results for N = ±2 
mod 7 given by (|83|) in (a) and by ((84l) in (b). Solid black lines represent the results 
for N = ±3 mod 7 determined by ([86| in (a) and by ([87j) in (b). 



To compare these formulas with the results of numerical simulations it is more 
precise to use the integrated nearest-neighbour distribution f|T6l) . In Figs. [T2] and [T3l 
such comparison is performed for all cases considered. The agreement is quite good and 
the differences are of the same order as in Fig. [2l 
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a) 



b) 



Figure 12. Difference between the integrated nearest-neighbour distribution and the 
corresponding theoretical prediction for a = f/5. (a) For non-symmetric ensembles, 
(b) For symmetric matrices. In each graph pictures differ by the matrix dimensions. 
From bottom to top TV = 801 {N = 1 mod 5) and 802 {N = 2 mod 5). 






a) 



b) 



Figure 13. The same as in Fig. [T2]but for a = 1/7. From bottom to top N = 799 
{N = 1 mod 7), 800 (A^ ee 2 mod 7) , and 801 (iV ee 3 mod 7). 



7. Summary 

A unitary random matrix ensemble 

Mkp = e' " _ ^2ni{k-p+aN)/N^ (^^) 

corresponding to a quantization of a simple pseudo-integrable map ([1]) is considered in 
detail. These matrices are characterized by a rational parameter a = m/q, the matrix 
dimension N, and symmetry properties of random phases $fc (jlj). To get a well defined 
limit of the spectral statistics of these ensembles for large it is necessary to consider 
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increasing sequences of such that niN has a fixed residue modulo the denominator 
of a 

mN = k mod q (89) 

with the residue k = 0,1, . . . ,q — 1. For /c = all eigenvalues of the main matrix 
can be found analytically as in [T3]. For all other residues the spectral statistics of 
the considered ensembles is non-trivial and differs considerably from standard random 
matrix ensembles. The cases k = 1 and k = q — 1 have been investigated in [15] where 
it was shown that for these k the nearest-neighbour distribution has the following form 

I a- '^e~'^^ for non — symmetric ensemble 
^i3/2-ig-(}s/2 symmetric ensemble 



In the present paper a kind of transfer operator method is developed to calculate the 
spectral statistics of the same matrix for all values of k. It is demonstrated that the 
nearest-neighbour distribution equals the product of a finite polynomial in s for non- 
symmetric matrices and in ^/s for symmetric matrices times the same exponential factor 
as in (|90|) 

k(q-k) 

anS^'e-'" for non — symmetric ensemble 

P(^) = { ikiT-k)-!) ■ (91) 

a„/2'S"''^e~'''^^^ for symmetric ensemble 

n=l 

Statistical properties of sequences with residue k and q — k are the same. The values 
of coefficients a„ can be calculated by finding the largest eigenvalue and corresponding 
left and right eigenvectors of the transfer matrix. 

For a = 1/5 and a = 1/7 and all possible residues the explicit form of these 
coefficients have been calculated. Numerical simulations in these cases are in a good 
agreement with obtained formulas. Other correlation functions can also be expressed 
from the same quantities. 

It appears that the considered ensembles of random matrices permit different 
generalizations which will be discussed elsewhere [19] . 
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Appendix A. 



The purpose of the Appendix is to give the proofs of certain formulas used in the text. 
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Let for all n = 1, . . . , 

m=l 

Solutions bjn of these equations can be expressed in terms of Cauchy determinants and 

Ylix^ - Vn) 
I I \-^m -^s ) 

A simple way to check it is to consider the function 

fn{x) = = . (A.3) 

s s 

Asymptotically fn{x) —>■ so the integral over a big contour encircling all poles equals 
1. Rewriting this integral as the sum poles gives 

YliXm-Vr) 

i = — ^Tnv ^ 

m \Xm yn) I I l^-^m Xg) 

which proves (IA.2p . 
Denote 

g{x) = . (A.5) 

x[[{x - Afc) 

k 

This function decreases as 1/x for large x and has poles at x = and x = A^. Integrating 
it over a contour encircling all poles one gets 

1 = e^™ + J2 ^^^j=F (A.6) 

m=l J_ J_ (Am — Afc) 

from which it follows that |-BmP defined by ( HTj) obey automatically the normalization 
condition (l33l) . 
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